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Abstract 



Starting from the Schrodinger functional, we give a non-perturbative definition of the 
running coupling constant in QCD. The spatial boundary conditions for the quark 
fields are chosen such that the massless Dirac operator in the classical background 
field has a large smallest eigenvalue. At one-loop order of perturbation theory, we 
determine the matching coefficient to the MS-scheme and discuss the quark mass effects 
in the /3-function. To this order, we also compute the Symanzik improvement coefficient 
necessary to remove the 0(a) lattice artefacts originating from the boundaries. For 
reasonable lattice resolutions and the standard Wilson action, lattice artefacts are found 
to be only weakly dependent on the lattice spacing a, while they vanish quickly with 
the improved action of Sheikholeslami and Wohlert. 
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1 Introduction 

The strong coupling constant a s can be extracted from experimental data, e.g. by com- 
paring jet production rates at a collider experiment. On the other hand, several authors 
have emphasized that a theoretical determination of the strong coupling constant is pos- 
sible by making use of lattice gauge theory [1-5]. This would provide a quantitative 
test of QCD which is believed to be the fundamental theory of both, hadronic physics 
and jet physics at high energies. 

The basic strategy, as proposed by Liischer, Weisz and Wolff [1] is easily de- 
scribed [6]. First, one fixes the free parameters in the QCD action by taking a cor- 
responding number of low energy observables as experimental input. Then one has to 
find a suitable non-perturbative definition of a running coupling and trace its evolution 
from low to high energies, using Monte Carlo simulations. In the high energy regime, 
perturbation theory can be used to convert to other schemes such as the modified min- 
imal subtraction scheme (MS) of dimensional regularization. 

Since the above strategy applies to any asymptotically free field theory, it has first 
been tested for simpler models, namely the 2- dimensional non- linear 0(3) model [1], and 
the pure SU(2) and SU(3) gauge theories [2-4]. One uses non-perturbatively defined 
couplings which run with L, the linear extension of the space-time volume. Then, using 
a recursive finite size scaling technique, the authors of refs. [1-4] were able to avoid 
the potential problem of having very different length scales fit on a single lattice. As 
a result, the running coupling was obtained deeply in the perturbative region, with all 
errors under control. The final perturbative matching to the coupling constant of the 
MS scheme is done at the one-loop level in SU(3), and in the case of SU(2) even the 
two-loop coefficient is known [7]. 

Although the choice of the non-perturbative running coupling is not a question 
of principle, the practical feasibility of the project very much depends on it. Indeed, 
the great success in the pure gauge theories is largely due to a clever definition of 
the running coupling constant as the system's response to an external color electric 
field. The theoretical framework for such a definition goes under the name Schrodinger 
functional. By this one means the Euclidean path integral on a space-time manifold with 
boundaries at Euclidean times xq = and xq = L, at which the values of the quantum 
fields are prescribed. The classical "path" then corresponds to a minimal action field 
configuration, the background field, which interpolates between the boundary values. 

Before the Schrodinger functional can be used in practice, some theoretical con- 
siderations have to be made. In particular, the presence of the boundaries may lead 
to new singularities, which are not taken into account by the usual renormalizations 
of the bare parameters. Moreover, the lack of translation invariance on the space time 
manifold prevents the use of the standard proofs of perturbative renormalizability. 

In the pure SU(N) gauge theory, these questions have been treated in ref. [8]. For 
a detailed account of the Schrodinger functional in QCD the reader should consult 
refs. [9,10]. Here we merely state the main result. Using power counting arguments, 
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one expects [11] that the Schrodinger functional in QCD is renormalized by the usual 
QCD renormalizations of the coupling constant, the quark masses and, in addition, a 
multiplicative renormalization of the quark boundary fields. In particular, if the latter 
are taken to vanish, no additional divergence is introduced through the presence of the 
boundaries. This naive expectation has been confirmed at one-loop order of perturbation 
theory [10]. 

In this paper, we adapt the definition of the SU(3) running coupling constant to 
QCD and establish its perturbative relation to the MS scheme at one-loop order of 
perturbation theory. The computation is straightforward in principle, but complicated 
through the presence of the quark masses. In particular, the dependence of the conver- 
sion coefficient and the /3-function on the quark masses cannot be computed analytically 
and has to be extracted numerically from lattice perturbation theory. In view of the 
Monte Carlo simulations to be carried out later, we also determined the lattice artefacts 
in the observables. 

The paper is organized as follows. In section 2, the definition of the Schrodinger 
functional is recalled both, in the continuum and on the lattice. The running coupling 
and the associated /3-function are defined. In section 3, the conversion to the MS scheme 
is done with and without the Sheikholeslami-Wohlert term in the action. There, we 
also discuss the /3-function and the threshold effects in the perturbative running of the 
coupling. Lattice artefacts are discussed in section 4 and we finally summarize our main 
results. 

2 The Schrodinger functional 
2.1 Classical continuum action 

A rigorous definition of the QCD Schrodinger functional exists on a space-time lattice 
and will be presented in section 2.6. On a more formal level, the Schrodinger functional 
can also be defined in the continuum [8,10]. It is given as the euclidean path integral 
with the action 



S[A,^] = - 




Here, go denotes the bare coupling constant, F^ v is the field tensor associated with the 
SU(3) gauge field A„, 

Fyi/ = d^Ajy — d^A^ -\- [Ay_, j4.j,], (2*2) 

and = + A^ denotes the covariant derivative on the quark fields. For simplicity, 
we assume rif degenerate quark flavors of bare mass mo, the generalization to the non- 
degenerate case being trivial. The 7-matrices are hermitian and satisfy 

{7^74 = 2^ W = O,-.. ,3, (2.3) 
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and the fifth 7-matrix is 75 = 70717273- Space-time is taken to be a cylinder of linear 
extension L. In the time direction, we impose the boundary conditions, 

Ak\ Xo =o = C k , A k \ X0=L = C' k , A: = 1,2,3 (2.4) 

and 

p+^L=o = 0, P-^\ X0=L = 0, 

(2-5) 

V>P-L 0=0 = o, ^+L 0=L = o, 

with the projectors P± = |(1 ± 70). In the spatial directions, periodic boundary con- 
ditions are imposed on the gauge fields, while the quark fields are taken to be periodic 
up to a phase factor, 

tp(x + Lk) = e i6 tp(x), tp(x + Lk) = e- i6 tp(x). (2.6) 

Here, the spatial index k takes the values 1,2,3, and k denotes the unit vector in 
direction k. For reasons to be explained in section 2.7, we will later consider the values 
and x/5 for the parameter 9. 



2.2 The background field 



The Schrodinger functional Z is considered a functional of the boundary gauge field C 
and C. In the following we restrict attention to the abelian boundary fields which have 
been introduced in ref. [3], 



C, = - 





with 



<t>i 

<t>2 
<t>3 



<t>2 



-<t>3 



2 )+ 3' 



<t>3 = -<t>2 



k — 1 ^ 2 , 3 5 



4x 

~3~' 
2x 

~3~' 
2x 

T' 



(2.7) 



(2.8) 



This defines a 2-parameter family of boundary gauge fields. A solution to the field 
equations with these boundary fields is given by 

#0 = 0, B k =[x C' k + (L-x )C k ]/L, As = 1,2,3. (2.9) 

Moreover, for given boundary fields, C and C, the field B represents the unique absolute 
minimum of the action. Any other field with the same action and boundary values is 
thus gauge equivalent to B, which will be referred to as the background field in the 
following. The associated field tensor G^ v has the non-vanishing components 

G ok = d B k = (C' k - C k )/L, k = 1, 2, 3, (2.10) 

which constitute the color-electric background field mentioned in the introduction. 
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2.3 The effective action 



The uniqueness of the induced background field allows to unambiguously define the 
effective action of the Schrodinger functional as a functional of B, viz 

T[B] = -In Z[C',C}. (2.11) 

In ref. [10], the saddle point expansion of the Schrodinger functional has been carried 
out to one loop-order, using dimensional regularization. At this order one needs the 
fluctuation operators of the ghost, gluon and quark fields. For the precise definition of 
the pure gauge theory operators, Ao and Ai, we refer to ref. [8]. 

Concerning the quark field fluctuation operator, it has been noted in ref. [9] that 
the Dirac operator 

V = D+m , D^ = 8^ + B^, (2.12) 

and its adjoint, , allow for the definition of two distinct operators 

A 2 = VV\ A' 2 = V ] V. (2.13) 

The operator A 2 is defined on spinors ?p(x) which satisfy the boundary conditions (2.5). 
Furthermore, its eigenfunctions satisfy the modified Neumann conditions 

(D - m )P_iP(x)\ X0=0 = 0, (D + m )P + iP(x)\ X0=L = 0. (2.14) 

On the other hand, the operator A 2 acts on fields which satisfy the same boundary 
conditions as t\) (2.5) and thus lives on a different space of functions. However, it 
follows from the analysis of ref. [9] that both operators, A 2 and A 2 , have exactly the 
same spectrum. For quantities which only refer to their spectrum we therefore need not 
distinguish between the two. 

To write down the one-loop effective action, we make use of the ^-functions, defined 
through 

C(s\Ai) = TrAr s , i = 0,1,2. (2.15) 

They extend to meromorphic functions in the whole complex plane, and one may show 
that their derivatives at s = are well-defined. One then obtains 

MsV 7 (2.16) 
- ^(01 Ax) + C'(0|A o ) + ^C'(0|A 2 ) + O(g^) , 

where g^r denotes the renormalized coupling constant in the modified minimal scheme 
(MS), and To is the classical action of the induced background field, 

ro[5] = -\ [ L d 4 xt T {G^G^}. (2.17) 
l Jo 
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Note that the operators A; have dimension L~ 2 . Therefore, the derivative of the (- 
functions contains a logarithmic dependence on L, such that a logarithm of fiL is ob- 
tained in eq. (2.16). Furthermore it turns out that the field dependent part of the 
one-loop effective action is proportional to IV 



2.4 Definition of the running coupling 

As the effective action depends - apart from the quark masses - only on a single scale, 
L, eq. (2.16) suggests to define a renormalized coupling constant, g(L), through [8] 



dv 

dr) 



9 2 {L) ' dr, 



= 12x. (2.18) 

?7=0 



Note that the derivative with respect to the parameter r, eliminates any divergent con- 
tributions to the effective action which do not depend on the background field. 
Using the notation 

we are here interested in the perturbative relation 

a m = a + Cl a 2 + 0(a 3 ) . (2.20) 

The coefficients in this expansion are functions of the parameter 

z = m(L)L, (2.21) 

where m is a suitably defined running quark mass. The form of these functions depends 
on the definition of m. However, for the computation of the one-loop coefficient c\ it 
is sufficient to define m at tree level. In the following we adopt the convention that, to 
this order, m coincides with the bare quark mass. 

We split ci into its pure gauge theory value, c^o, and the quark contribution, c\ t \, 

viz 

ci = ci(n/, z) = ci, + n f c ltl (z), (2.22) 
where c^o has been computed in ref. [3], 

c 1]0 = 1.25563(4). (2.23) 



different from zero, viz 



A further renormalized quantity is obtained if the parameter v in eq. (2.18) is kept 

{j2-™\ (2.24) 



dr] 



= 12?H 

?7=0 
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It is not difficult to see that v does not depend on v. Note also that v vanishes at tree 
level, because To is independent of v. Its perturbative expansion reads 

v = t>i(n/, z) + 0(a), t>i(n/, z) = v 1>0 + rifV ltl (z), (2.25) 

with the pure gauge theory contribution [3] 

«i,o = 0.0694603(1) . (2.26) 

The coefficients c\ and v\ can be calculated in the continuum through evaluation of 
eq. (2.16). One starts from the expression [8] 

-C'(0|A) = lim{( 7E + m*)-a (A) - £ e"^ — lnA n }, (2.27) 

where A stands for any of the fluctuation operators. The Seeley coefficient ao is pro- 
portional to To in all three cases, and A n are the eigenvalues of A, in ascending order. 

The calculation of the quark contribution c\^z) essentially amounts to the deter- 
mination of the eigenvalues of A2 and their derivatives with respect to tj, up to a certain 
cutoff in the level n. One then evaluates the bracket in eq. (2.27) for a range of ^-values 
and extrapolates to 8 = 0, taking into account that the bracket has an asymptotic 
expansion in powers of 8 1 ! 2 [8]. 

We employed a variational method with a plane wave basis to compute the eigen- 
values. As a check on the precision, we also considered the eigenvalue equation for A2, 
which has a general solution in terms of hypergeometric functions. The boundary condi- 
tions eqs. (2. 5), (2. 14) then lead to a system of linear equations which can be numerically 
solved for the eigenvalues. 

For z = 0, we were thus able to compute ci^ with an estimated numerical precision 
of 3-4 significant digits. As z increases, cancellations in the sum over eigenvalues become 
stronger, resulting in a loss of precision. Since this computation was mainly intended 
to be a check on the more precise lattice methods, we do not quote the results here and 
refer the reader to section 3 instead. 



2.5 The /3-function 

The Callan Symanzik /3-function is defined through 

^) = -i||, (2.28) 

where the derivative is taken at fixed bare parameters or, equivalently, keeping the 
renormalized parameters, m(Lo) and g(Lo), fixed at some normalization scale Lq. The 
/3-function has a perturbation expansion 

OO 

m S -° ~f E h nf n , (2-29) 

n=0 



6 



with coefficients that are quark mass dependent, 

b n = b n (n f ,z), n= 0,1,2,... (2.30) 

In particular, 60 and b\ coincide with the universal coefficients only when the quark 
mass is set to zero, 

6o(n/ ' 0) = (4^ (11 "I n/) ' (2l31) 
1 38 

6l(n/ ' 0) = (4xF (l ° 2 "T n/) - (2l32) 

To obtain the mass dependence of the coefficients b n , one may relate g perturbatively to 
any mass independent renormalized coupling constant. The coupling in the MS scheme 
has this property. Writing 

b (nf, z) = 6 ,o + nfb 0tl (z), (2.33) 
we obtain 6o,i(z) from ci^z) through 

Vi(^) = "2^2 - ^ zc 'i.i( z )' c 'i.i( z ) = Tz Cl ' l{z) - (2-34) 

The central observable in a non-perturbative computation of the evolution of the 
coupling is an integrated version of the /3-function, the step scaling function [1]. It is 
defined as follows. Starting with a value u = g 2 (L) for the coupling at length scale L, 
the step scaling function, a, is 

a(s,u,z) = f(sL) . (2.35) 

Finally, we note that a similar scaling function can be defined for the running quark 
mass. 

2.6 Lattice formulation 

2.6.1 The lattice action 

In Wilson's lattice QCD, the path integral representation of the Schrodinger functional 
reads [9], 

Z = J~D[iP]D$]D[U]e- S . (2.36) 
with the lattice action S = S g + Sf, given by 

S 9lU] = 4 E™(p)tr{l -[%)}, (2.37) 
9o p 

S f [U,$,il>] = a^^D + mo)^- (2.38) 
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The gauge field action S g is a sum over all oriented plaquettes p on the lattice, with the 
weight factors w(p), and the parallel transporters U(p) around p. The weights w(p) are 
1 for plaquettes in the interior and 

I |c s if p is a spatial plaquette at xq = or xq = L, 
W \P) = \ 

yc t if p is timelike and attached to a boundary plane. ^ 39) 

The choice c s = c t = 1 corresponds to the standard Wilson plaquette action. How- 
ever, these parameters can be tuned in order to reduce the lattice artefacts, as will be 
discussed in detail in section 4. 

The Dirac operator in the quark action (2.38) is specified by 

D = \ E{7m(V; + V„) - aV^V^} + c sw j £ (2.40) 

fi=0 i_i,i'=0 

with the forward and backward covariant derivatives 

V^z) = hu(x,n)iP(x + ajl)-iP(x)}, (2.41) 

V*V(z) = ^[ip(x)-U(x-afi,nyip(x-afi)]. (2.42) 

The link variable U(x, /x) is the usual parallel transporter from point x + afi to point x, 
where fi denotes the unit vector in /x-direction. 

The last term in eq. (2.40) has been introduced by Sheikholeslami and Wohlert [12], 
in order to cancel the leading cutoff effects of the standard Wilson quark action. It 
contains the lattice definition P^ v of the field tensor, 

P^A X ) = ^| U(x,fi)U(x + afi,v)U(x + av, fi)^U(x,v)^ 

+ U(x, v)U(x + av — afi, fi)^U(x — afi, v)^U(x — afi, /x) 

+ U(x — afi, [i)^U(x — afi — av, v)^U(x — afi — av, fi)U(x — av, v) 

+ . +1 ( 2 - 43 ) 

+ U(x — av, vyU(x — av, fi)U(x -\- afi — av, v)U(x, /x) T > 

- (/x < — > v), 
and our convention for reads 

<V< = ^[7/x,7k]- ( 2 - 44 ) 

The standard Wilson quark action is recovered for c sw = 0, and the choice c sw = 1 will 
be referred to as Sheikholeslami- Wohlert action. 
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2.6.2 Boundary conditions and the background field 
The boundary conditions for the lattice gauge fields are 

U(x,k)\ X0=0 = exp(aC fe ), U(x,k)\ Xo=L = exp(aC£), (2.45) 

for k = 1,2, 3, and with the abelian boundary fields C and C as given in eq. (2.7). All 
other boundary conditions are as in the continuum (cf. section 2.1). 

The boundary conditions (2.45) lead to a unique (up to gauge transformations) 
minimal action configuration V, the lattice background field. It can be expressed in 
terms of B (2.9), 

V(x,n) = exv{aB^(x)}, (2.46) 
and its field tensor (2.43) evaluates to 

P ^\u=v= ^2 sinh « 2 GW, ( 2 - 47 ) 

with G^ v given in eq. (2.10). 

2.6.3 The effective action 

The effective action T = — In Z has an asymptotic expansion in the bare coupling 
constant, 

r = go 2 T + I\ + O(g 2 ) (2.48) 

with the lowest order term To = {fl , o^g[^]} go - - ^he nex t order term is, for c s = c t = 1, 
given by 

Ti = ^ In det A x - In det A - ^ In det A 2 . (2.49) 

Here, the operators A;, i = 0, 1, 2, are the lattice approximants of the continuum oper- 
ators introduced above. Again, we refer to refs. [8,3] for the definition of Ao and Ai. 
The operators A 2 and A' 2 are related to the lattice Dirac operator, 

A 2 = [(D + m ) l5 ] 2 , A' 2 = [ 7B (D + m )] 2 . (2.50) 

Since D (2.40) acts on lattice spinors ^(x) which satisfy the boundary conditions (2.5), 
it follows that the eigenfunctions of ~f5(D + mo) (and thus of A 2 ) satisfy a lattice version 
of the modified Neumann conditions (2.14). Furthermore, from eq. (2.50) it is obvious 
that both lattice operators, A 2 and A 2 , have the same eigenvalues. 
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2.6-4 Definition of the running coupling 



With these preliminaries, we define the running coupling constant g(L) and its relative 
v(L) through 



dv 

dn 



k 



T) = 



with 



dT 



dr) 



T) = 



12(i/a) 2 [sin( 7 ) + sin(2 7 )], 7 = -x(a/i) 2 . 



(2.51) 



(2.52) 



Here, the normalization constant k ensures that the renormalized coupling is exactly 
equal to go at lowest order in the perturbation expansion. 



2.7 Spatial boundary conditions and the Dirac operator 

It remains to justify the choices and x/5 for the angle 9 in the spatial boundary con- 
ditions for the quark fields (2.5). Of course, setting 9 = is a natural and aesthetically 
pleasing choice, which corresponds to periodic boundary conditions. 

Beyond aesthetical criteria there is an additional, more technical one. Having in 
mind numerical simulations of QCD, it is worthwhile to recall that the speed of the 
known algorithms depends crucially on the condition number, i.e. the ratio between the 
highest and the lowest eigenvalue of the squared fermion matrix A2 (2.50). Thus it is 
desirable to have a lowest eigenvalue of A2, which is not too small. 





9 = 








9 = x/5 






n 




n c 


d 


n 




n c 


d 


1 


2.132449 


2 


2 


1 


4.693976 


2 


2 


2 


4.804360 


2 


2 


2 


4.881719 


1 


2 


3 


7.599922 


3 


2 


3 


8.384625 


2 


2 


4 


9.732686 


1 


2 


4 


13.109607 


1 


2 


5 


12.132366 


3 


2 


5 


13.834839 


3 


2 


6 


20.221614 


1 


2 


6 


19.717052 


3 


2 


7 


22.937850 


2 


2 


7 


26.100020 


2 


2 


8 


23.755510 


2 


2 


8 


26.996453 


2 


2 


9 


27.184931 


1 


6 


9 


27.826014 


1 


6 


10 


28.451394 


3 


6 


10 


27.846538 


3 


6 



Table 1: The lowest eigenvalues (in units of L~ 2 ) of the continuum operator A 2 for vanishing mass 
and two choices of 9. Since A£ is diagonal in color space, each eigenvalue is associated with a color 
component n c . We also give the degeneracy d for one quark flavor. 

As a guiding principle, we have investigated the eigenvalues of the corresponding 
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continuum operator 1 A£ (2.13). In ref.[9], it has been emphasized that the boundary 
conditions (2.5) introduce a gap, i.e. a non-vanishing lowest eigenvalue which is given 
by 7r 2 /4£ 2 = 2.467. ../L 2 , if the phase 9, the quark mass and the background field are 
taken to vanish. 



9 = 




L = 


6a 


L = 


12a 


L = 


24a 






n 


^sw — 1 


c sw — 


^sw — 1 


c sw — 


^sw — 1 


c sw — 


n c 


d 


1 


2.5338 


3.0378 


2.3147 


2.5504 


2.2190 


2.3330 


2 


2 


2 


5.7261 


5.0256 


5.2428 


4.9024 


5.0181 


4.8503 


2 


2 


3 


9.4614 


10.4139 


8.4086 


8.8526 


7.9758 


8.1907 


3 


2 


4 


12.5904 


14.8258 


10.9562 


11.9796 


10.2964 


10.7871 


1 


2 


5 


14.4516 


13.3254 


13.1233 


12.5815 


12.5924 


12.3256 


3 


2 


6 


23.6179 


20.7726 


21.4968 


20.1181 


20.7701 


20.0870 


1 


2 


7 


25.3644 


26.1068 


24.3912 


25.1434 


23.7033 


24.0727 


2 


2 


8 


27.6740 


26.9095 


25.9626 


25.1982 


24.8980 


24.5223 


2 


2 


9 


32.0603 


33.9463 


29.1583 


29.7601 


28.0767 


28.2945 


1 


6 


10 


33.4272 


34.4403 


30.7130 


31.0696 


29.5507 


29.6932 


3 


6 


9 = x/5 


1 


5.7398 


6.4868 


5.1553 


5.5042 


4.9100 


5.0788 


2 


2 


2 


6.0947 


7.6777 


5.4075 


6.1333 


5.1253 


5.4732 


1 


2 


3 


9.9756 


9.0438 


9.0987 


8.6484 


8.7240 


8.5022 


2 


2 


4 


15.2653 


12.9957 


14.0326 


12.9231 


13.5401 


12.9897 


1 


2 


5 


17.6812 


18.9750 


15.4669 


16.0670 


14.5849 


14.8749 


3 


2 


6 


23.6820 


22.2339 


21.2686 


20.5772 


20.4024 


20.0625 


3 


2 


7 


29.9843 


30.9829 


28.0259 


28.8318 


27.0343 


27.4287 


2 


2 


8 


32.6882 


31.6622 


29.8174 


28.9978 


28.3772 


27.9760 


2 


2 


9 


31.9457 


33.8957 


29.7929 


30.4575 


28.7990 


29.0545 


1 


6 


10 


33.8957 


34.8946 


30.2866 


30.6210 


28.9484 


29.0768 


3 
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Table 2: The lowest eigenvalues (in units of L 2 ) of the lattice operator A2, for vanishing bare mass 
mo. The eigenvalues are ordered according to their continuum limits (cf. table 1). 

This picture remains valid in the presence of the abelian background field (2.9). 
At 9 = and for vanishing quark mass, the lowest eigenvalue is slightly decreased and 
approximately given by 2.132/i 2 . Furthermore, it is possible to increase the lowest 
eigenvalue substantially by varying the angle 9. We observed a maximal gap around 
9 = x/5, and decided to consider this value besides 9 = 0. The lowest eigenvalues of 
for these two choices are shown in table 1. 

x In this subsection we add a superscript c in order to distinguish the continuum operator from its 
counterpart on the lattice. 
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On a lattice with finite spacing a, the gap is somewhat larger. For illustration, 
we have collected the lowest eigenvalues of the lattice operator A2 on three lattices 
of different size L/a in table 2. We observe that almost all low lying eigenvalues are 
approached from above by their lattice approximants, at a rate which is roughly given 
by a/L. In this respect, there is no difference whether the Sheikholeslami-Wohlert term 
is included in the action or not. We notice, however, that the approach to the continuum 
eigenvalues seems to be slightly more uniform for the Sheikholeslami-Wohlert action. 

It is not clear how relevant these considerations are for a realistic simulation. Fluc- 
tuating gauge fields can lead to small eigenvalues of the lattice operator A2 and finally 
the typical spectrum has to be determined by numerical experiment. However, we be- 
lieve that in sufficiently small volume, the dominant contributions to the path integral 
come from gauge field configurations which are close to the classical background field. 
Thus we suggest to use 9 = x/5 in practical applications. 



3 One loop relations 

The one-loop matching between a and a^jg (2.20) is done in two steps. First, we 
calculate the relation of a to a renormalized lattice coupling ai at . Combination with 
the known one-loop relation between ai at and a^nr then yields the desired result. 



3.1 The basic calculation 

Insertion of the asymptotic expansion (2.48) into the definition of the running coupling, 
eq. (2.51), leads to the relation 

f = 9o+Pi9o + 0(g 6 ), (3.1) 
with pi being explicitly given by 

1 dTt 



Pi 



k drj 



(3.2) 

?7=f=0 



The coefficient p\ depends on the number of flavors, the bare quark mass in lattice units, 
mod, and the lattice size I = L/a. For later convenience, this dependence is written in 
the form 

Pi = Pi(nf, z, I) = Pi,o(0 + n f p ltl (z, I), (3.3) 

where we have set 

z = JnL, Tn = — ln(l + mod). (3-4) 

To this order, we thus identify m with the so-called pole mass, which is in one-to-one 
correspondence with the bare quark mass and coincides with the latter to leading order 
in the small a expansion. 
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When the lattice spacing becomes small, the field dependent part of the one-loop 
effective action, Ti, is expected to diverge logarithmically. As usual, this divergence 
is absorbed in the renormalization of the coupling constant. We thus eliminate go m 
favor of a renormalized coupling g lat (at the renormalization scale /x), defined through 
minimal subtraction of logarithms [13], i.e. 

9o = 9lat + z i( n f, a H)9lat + O(ffL), (3.5) 

with 

zi(nf, a/x) = zi ]0 (a/i) + ra/Zi^a/x) = 26 (n/, 0) ln(a/x). (3-6) 
Eq. (3.1) then becomes a relation between renormalized coupling constants, 

t = dlt + (Pi + zi)9L + O(fff ot ), (3.7) 

and we expect the coefficient in the r.h.s. of eq. (3.7) to have a well-defined continuum 
limit. Indeed, for rif = 0, this has been demonstrated in refs. [8,3]. 

Therefore, we may restrict attention to the analysis of the quark field contribution, 
which, for c s = c t = 1, is given by 

d 

Pi,i(z,/) = (2kn f )- 1 — lndetA^^. (3.8) 

As we are considering spatially constant abelian background fields, eq. (2.9), A2 is 
diagonal with respect to its color and spatial momentum dependence. Its determinant 
factorizes accordingly and the problem is reduced to the evaluation of the determinant 
of a 4(1 — 1) X 4(1 — 1) matrix in each subspace of fixed spatial momentum and color. 
In these subspaces, A2 acts as a second order difference operator in the Euclidean time, 
with coefficients that are 4x4 matrices in Dirac space. Following appendix C of ref. [8], 
the determinant of such an operator can be efficiently calculated by solving a recursion 
relation. Differentiation with respect to tj then leads to a recursion for the quantity of 
interest, i.e. the 7/-derivative of the determinant. 

Alternatively, one may use the fact that A2 is the square of a first order opera- 
tor (2.50). In the subspaces of fixed momentum and color, the determinant of ^^(D+mo) 
can be computed by solving a first order recursion relation. We include a more detailed 
discussion of this method in the appendix. 

We did the computation using either method and obtained pi^(z, I) for lattice sizes 
up to / = 64 and in "REAL*16" precision. For fixed value of z = mL, we expect pi t i to 
be given by an asymptotic series of the form [8] 

Pi.i ~ r + s ln/+ (ri + siln/)//+ 0(ln/// 2 ). (3.9) 

The first few coefficients in eq. (3.9) can be extracted efficiently by first cancelling higher 
order terms in l/l through numerical differentiation and then checking for stability as 
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I is increased [14]. Due to the large range of I and the good numerical precision, this 
technique gives very accurate results for fixed values of z. In this way we determined 
the coefficient so of the logarithm in eq. (3.9), 

5 = -0.00844344(2). (3.10) 

This value has been obtained for z = 0, 9 = x/5 and c sw = 1. Furthermore, within the 
numerical accuracy, we found this result to be independent of all three parameters. 

Recalling eqs.(3.6) and (3.7), we see that so coincides with the expected result 
s = 26 ,i(0) = -l/(12x 2 ) = -0.0084434319.... The continuum limit of the quark 
contribution to pi + z\ thus exists and is given by 

lim[pi i(z, L/a) + z 1 ,i(a/x)] = r (z) + s ln(/xi). (3.11) 

We will set /x = i _1 in the following and assume the exact value so = — l/12x 2 in the 
numerical analysis of the series (3.9). Since, in this section, we are interested in the 
universal relations between renormalized couplings in the continuum, we will concentrate 
on the conversion coefficient ro(z). The higher order terms in eq. (3.9) represent lattice 
artefacts and will be discussed separately in section 4. 

3.2 Matching coefficient for massless quarks 

The quark contribution to the matching coefficient between a and cti at = flf a i/47r is 
4xro(z). For massless quarks, tq can be read off from table 3. Note that ai at depends 
implicitly on the action, i.e. on the coefficient c sw . 








'o(O) 


x/5 


1 


-0.034664940(4) 


tt/5 





-0.009868186(4) 





1 


-0.03328359(1) 








-0.00848683(1) 



Table 3: The first term in the asymptotic expansion (3.9) of 

In order to obtain the one-loop relation between a and a^jg, we also need the 
relation of ai at to a^jg, 

a iat = «Ms + d i( n /) a Ms + °( a Ms)- ( 3 - 12 ) 

For the case of the Wilson action, the coefficient d\(rif) = c^o + 1S known in 

the literature [15,16]. In particular, the quark field contribution d\ t \ has first been 
calculated by Weisz [15], who gives the result in terms of the lattice integral P3, 

d 1A = 4xP 3 , (3.13) 
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which was determined numerically with an accuracy of 3-4 significant decimal places. 
A more precise value for P3 can be extracted from ref. [16], where a precision of 6 
significant digits was achieved. We have re-evaluated P3 and obtained 

P 3 = 0.0066960(1), (3.14) 

which is, within errors, in agreement with the aforementioned values. 

Taking T"o(0) for c sw = and both values of 9 from table 3, we finally obtain the 
matching coefficients ci^ for massless quarks. They are both small and positive, 

f 0.039863(2) for 9 = x/5, 
c M =-4xP 3 + ^o0 = 7 ' 3.15 

[0.022504(2) for 9 = 0. 

As already mentioned in section 2, these values are in agreement with the ones obtained 
from the direct expansion of the effective action in a^jg. 



3.3 Massive quarks 

Let us now turn to the mass dependence of a, which, at one-loop order, is described by 
ci t i(z), eq. (2.22). With our numerical methods we can compute ci^z) for fixed values 
of z, but these methods do not allow for a direct determination of the z-dependence. We 
therefore chose to represent the z-dependence of ci ] i(z) by a Chebyshev approximation 
(T n (x) = cos[raarccos(a:)]) 

ci,i(z)-c 1 , 1 (0) = z[^i n T n (z/5-l) ±10" 6 }, zG[0,10]. (3.16) 

The coefficients t n are listed in table 4. For most purposes, it is sufficient to use a 
truncated version of this Chebyshev approximation. For example, the first five terms 
describe ci ] i(z) — ci^O) with an accuracy of 5 • 10 _4 z. 

For large z, one needs a different approximation to ci i(z). In order to arrive 
at an appropriate form, we note that a heavy fermion is expected to decouple from 
any physical prediction in the limit z — ► 00 [17]. An example is the relation between 
physical (and thus quark mass dependent) couplings. The fermionic contribution to the 
perturbative matching coefficients therefore should vanish for large values of z, up to 
power corrections in 1/z. On the other hand, in unphysical definitions of the coupling, 
such as the MS scheme, heavy quarks do not decouple. Therefore, we expect ci ] i(z) 
to be logarithmically divergent in the large z limit. Furthermore, the above argument 
shows, that not only the logarithm, but also a subleading constant contribution to ci^ 
would be the same in the relation between a^jg and any other physical coupling. In 
particular, in the case of the coupling defined via the static quark potential, this constant 
contribution is absent. We are therefore led to approximate 

1 3 

ci,i(z) = -— ln(z)+ ^g n z- n ±5-10" 5 , z" 1 < 0.13 , (3.17) 
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n 












= 7T/5 




V — x/5 


= 





-0.05996947 


-0.05560781 


C\ C\C\ O O C\Oi V 

-0.0066985 


-0.0066991 


1 


0.04146952 


0.03840671 


0.0057316 


C\ C\C\Ci A A i"k O 

0.0084403 


2 


-0.01440583 


-0.01451262 


-0.0017853 


-0.0045988 


3 


0.00411117 


0.00535508 


0.0002527 


C\ A A A A A A V 

0.0020945 


4 


-0.00068390 


-0.00187650 


0.0001202 


-0.0007739 


5 


C\ f\f\f\~\ O AAO 

-0.00013008 


rk rk f\ f\ o s\ c\ o c\ 

0.00062089 


-0.0001110 


0.0002192 


6 


0.00016711 


-0.00019288 


0.0000473 


r\ AAAAO A1 

-0.0000391 


7 


-0.00007951 


0.00005478 


-0.0000112 





Q 
O 


U.UUUUZ04 i 


n nnnm oqq 
— u.uuuuizyy 






9 


-0.00000288 


0.00000168 






10 


-0.00000143 


0.00000059 






11 


0.00000107 


-0.00000049 







Table 4: Chebyshev coefficients of the fits eqs. (3. 16), (3. 20). 



with 

/(-0.10822, -0.0013,0) for 9 = x/5, 

(?i5 ?2, q3) = < (3.18) 

[ (-0.10718, -0.0334, 0.233) for 9 = 0. 

The coefficients (i = 1, 2, 3), have been obtained by fitting eq. (3.17) to a number of 
data points ci^z;) with 6 < z; < 25. To assess the quality of the fit, we varied the 
number of terms in eq. (3.17), and the number of points z; used for the fit. We conclude 
that ci ] i(z) — (^1(0) is represented with a precision of 5 • 10~ 5 for z _1 < 0.13. 

3.4 Result for v 

The quantity v defined in eq. (2.51) is obtained in non-perturbative Monte Carlo sim- 
ulations without any extra computational effort. Therefore it represents a useful check 
of the applicability domain of perturbation theory and the approach to the continuum 
limit [3]. 

In perturbation theory, we determined the contribution of massless quarks to 
eq. (2.25), 

f 0.0245370(1) for0 = 7r/5, 
v hi = < K ' 3.19 

[0.013554(1) for = 0. 

Again the quark mass dependence can be accurately described using a Chebyshev fit, 

. . f «i,i(0) + z El =0 ^nT n (z/5 - 1) ± 5 • 10" 5 ] for z £ [0, 10], 
V\ i(z) = < 

[O ± 5 • 10- 5 for z > 8 ( 3 20 j 

with the coefficients w n as given in table 4. 
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3.5 A-parameters 

The one-loop coefficients between different renormalized couplings translate to ratios of 
the corresponding A-parameters. More precisely, if the one-loop relation between the 
couplings g a and g b reads 

gl = gl + c ab gt + 0(g 6 b ), (3.21) 

the ratio between the A-parameters is given by 

A a /A b = e c ^ 2b °. (3.22) 

Here both, 60 and c ab refer to the coefficients for vanishing quark masses. The reason 
is that A-parameters refer to the asymptotic high energy regime where all mass effects 
are negligible. 

In the pure SU(3) gauge theory, the ratio between the A-parameters associated 
with the Schrodinger functional and the MS scheme is 

A/A M g = 0.48811(1). (3.23) 

If the quark fields are included, this value is lowered. For instance, for 9 = x/5 and 
rif = 1 (rif = 3), the approximate value is 0.455 (0.383). Similar values are obtained 
with 9 = 0, and, recalling the relation between the A-parameters of the MS and MS 
schemes, 

A ms/Ams = ( e7E l^f' 2 = 0.37647475..., (3.24) 

we conclude that A is, for three or four quark flavors, almost equal to Ams- 

As a byproduct of our computational strategy we also obtain the relation between 
the A-parameters Asw an d Aw , associated with the lattice couplings, gi at , of the 
Sheikholeslami-Wohlert [12] and the standard Wilson action, 

ln(A S w/Aw) = n f (r (0)| Csw=0 - r (0)| Csw=1 )/26 

= 0.024796754(2) n f /2b . (3.25) 

This relation holds for an arbitrary number of colors N with 60 = (UN — 2n^)/48x 2 . 
This is easy to understand when the matching of the couplings is done using the back- 
ground field technique adapted to lattice perturbation theory [18,15,19]. Considering 
the 2-point function of the background field at one-loop order, the only difference be- 
tween the two schemes are the quark loop diagrams with at least one quark-gluon vertex 
stemming from the Sheikholeslami-Wohlert term in eq. (2.40). Since this vertex is pro- 
portional to a group generator in the fundamental representation, the contraction of the 
group indices yields a iV-independent constant. 

As a further check of this statement, we performed the same computation using 
the abelian SU(2) background field of ref. [8] and obtained agreement within the nu- 
merical accuracy. Indeed, the number quoted in eq. (3.25) is the one obtained from this 
computation. 
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3.6 Threshold effects in the /3-function 



According to eq. (2.34), the mass dependence of the one-loop /3-function can be obtained 
from ci t i(z). The derivative zc[ x (z) which appears in 6o,i(z) is, for large values of z, 
well approximated by taking the derivative of the fit eq. (3.17). On the other hand, for 
small values of z we may differentiate the Chebychev fit, eq. (3.16). In this way, we 
obtained zc[ x (z) with an estimated precision of 10~ 4 in the whole range of z- values. 




The resulting function &o,i( 2 ) (2.33) is displayed in figure 1. One notices that 
the transition from an effectively massless quark to an approximately decoupled heavy 
quark is not very rapid. The reason is the following. 6o,i(z) — 6o,i(0) starts with a term 
proportional to z at small z. At the other end, the decoupling limit, lim^oo &o,i( 2 ) = 
0, is approached with a 1/z- correction. These correction terms that are odd in the 
quark mass are possible because the boundary conditions eq. (2.5) do not respect chiral 
symmetry. This entails a rather broad transition region, in contrast to the case of other 
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mass-dependent schemes, such as the MOM-scheme [20] or the coupling defined via the 
static quark potential. 

At first sight, the existence of the 1/z-term in the /3-function of a heavy quark 
seems to pose a problem for the full non-perturbative computation of the /3-function 
along the lines of ref. [3]. The reason is that quarks that are much heavier than 1/L 
contribute lattice artefacts that dominate over their physical effect if one is limited to 
lattice sizes of, say, L/a < 20. One would therefore like to omit quarks with z > z cut 
from the Monte Carlo simulations, with a realistic z cut taking values between, say, 2 
and 4. 

Our result allows to quantify the error that is made in the one-loop evolution of the 
coupling g, if a quark with z > z cut is omitted in the /3-function. We fix the value of g(Lo) 
for some value mLo = zo ^ 1 and want to know the coupling when mL = z < 1. (One 
may think, e.g. of L = (0.5 GeV)" 1 , m = 4.5 GeV and L = (20 GeV)" 1 ). Omitting a 
heavy quark in the one-loop /3-function corresponds to setting 6o,i(z) = for z > z cut . 
The ensuing error is the difference in the couplings obtained by running once with the 
full and once with the truncated /3-functions from scale Lq to L. It is convenient to 
look at the corresponding error of g~ 2 which is given as the integral over the tail in the 
/3-function between z cut and zo (cf. fig. 1), 

A(ff" 2 ) = i-{ci,i(zo) - ci,i(z cut ) + -L ln(z /z cut )}. (3.26) 

Being generous, we enlarge the error by taking the limit of large zq. This corresponds to 
taking the integral over the whole tail of 6o,i(z), which approximately yields A(g~ 2 ) = 
0.003 for z cut = 2 and both 9 = and 9 = x/5. In comparison, the present experimental 
error for gZ?L(Mz) from the LEP-experiments is about 0.05 [21], and the pure gauge 
coupling is known with precision A(g~ 2 ) = 0.02 [3] (Note that the error of g~ 2 is 
preserved under one-loop evolution of the coupling). 

If the real error is of the order of magnitude suggested by one-loop perturbation 
theory, one may conclude that heavy flavors with, say z > 2, can safely be omitted from 
a Monte Carlo simulation. On the other hand, perturbation theory suggests that quark 
flavors with z < 2 do not cause significant cutoff effects as we will see in section 4. 
Provided this picture is correct also beyond one-loop perturbation theory, the evolution 
of g can be computed through Monte Carlo simulations including the effect of a quark 
flavor of any mass and with all errors under control. 

4 Lattice artefacts 

The evaluation of finite lattice spacing effects in perturbation theory is important in sev- 
eral respects. First, we gain some insight about the cutoff effects to be expected in a full 
non-perturbative computation. Second, they are needed to determine the (Symanzik-) 
improvement coefficients and third, it has proven to be useful to define perturbatively 
improved observables [4], for which the one-loop cutoff effects are cancelled completely. 
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4.1 Symanzik's local effective Lagrangian 

Consider a general lattice action for the Schrodinger functional. We expect that observ- 
ables, such as g, converge to their continuum values with dominant corrections that are 
linear in the lattice spacing a (up to logarithmic corrections). This expectation is derived 
from the following reasoning. According to Symanzik [22] we can represent the lattice 
theory with finite spacing a by a continuum theory, with a local effective Lagrangian 
containing the terms in eq. (2.1) plus higher dimensional operators accompanied by 
explicit powers of a. As long as these are not forbidden by an exact symmetry of the 
lattice action, we have to expect that all composite operators, which can be formed from 
the basic fields, are present in the local effective Lagrangian. In QCD with fermionic 
fields, gauge invariant dimension five operators exist. Integrated over the volume, they 
can contribute 0(a) lattice artefacts [12]. In addition, in the Schrodinger functional, 
dimension four operators integrated over the surfaces at xq = and xq = L are a source 
of 0(a) terms. 

0(a) perturbative improvement eliminates these leading lattice artefacts order by 
order in perturbation theory. This is achieved by adding lattice representatives of the 
corresponding operators to the action one starts from. Their coefficients can then be 
tuned to cancel all 0(a) effects to a given order in the coupling. 

In the pure gauge theory, dimension five operators do not exist. Gauge invari- 
ance and axis permutation symmetry restrict the set of dimension four operators to 
k=i ^ T {FjkFjk} an d Sfe=i tr{-^ofc-Fofc}- Our action (section 2.6) contains these oper- 
ators with a strength tunable by the coefficients c s and c t , respectively. For abelian 
background fields, as considered here, the first operator vanishes identically at the sur- 
faces. Consequently, we do not obtain any information on c s , here. In order to discuss 
perturbative improvement, we expand the second coefficient as 

et = 1 + #\n f )gl + . . . , «£>(»,) = + nfci 1 ^ . (4.1) 

Here, the 0(1) term is fixed by requiring that the linear term is absent in the small a 
expansion of To and the pure gauge part, 

c| 1,0) = -0.08900(5), (4.2) 

is known from ref. [3] . 

Once we include the fermion fields, we have to consider two operators of dimen- 
sion five [12]. Sheikholeslami and Wohlert have shown that one of these operators is 
redundant, when the theory posesses full translational invariance [12]. It is easy to see 
that the same is true for the Schrodinger functional. Consequently we include only the 
remaining operator, with its strength determined by c sw . Setting c sw = 1 corresponds 
to tree-level improvement. In the Schrodinger functional, one also needs to consider 
possible dimension four operators - such as do^tptp) - at the boundaries. While we 
postpone a complete analysis of these terms to future work, we conjecture that they are 
included in our action with the correct weight for tree-level 0(a) improvement. Indeed, 
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this is indicated by the cutoff dependence of the observables g and v, as will be seen 
below. 

4.2 0(a) improvement 

Let us step back and discuss eq. (3.9) in more detail. Its structure is dictated by 
the aforementioned equivalence of the lattice theory with an effective continuum theory 
containing local interactions [22]. For example, the operators discussed above contribute 
artefacts that are linear in the lattice spacing. At one-loop order, the loop integration 
will generally introduce an additional logarithmic dependence on the lattice spacing. 

These arguments lead us to expect that the continuum limit is approached in the 
form eq. (3.9) with the coefficients tq . . . being functions of z. For a tree-level improved 
action, s\ should vanish. Indeed, setting c sw = 1, we find s\ = within our numerical 
accuracy of about 10~ 4 and for both values of 6. 

We are left with the linear term ri(z)/Z, which we determined numerically (for 
c sw = 1) to be 

{0.038282(2) for 9 = x/5, 
0.038282(2) for = 0, (4.3) 
0.0382820(1) in the SU(2)-theory with = 0, 

and 

ri(z)-ri(0) = 0.012 z ± 2 • 10" 5 for z < 10 . (4.4) 

There are no local operators that could generate terms like z 2 /I or higher orders in z. 
This argument explains the structure of our numerical result (4.4). On the other hand, 
the linear term can be rewritten z/l = amo + 0(/~ 2 ), corresponding to the operator 
motT{F fJil ,F f j il ,} in the local effective Lagrangian. In order to remove it in the improved 
action one just needs to perform a redefinition of the bare coupling go. F° r this reason, 
this term does not represent a genuine lattice artefact: its effect vanishes as soon as 
relations between renormalized quantities are considered. An example is provided by 
the step scaling function a. Note that a similar reparametrization applied to the bare 
mass is expected to be necessary, when one considers the continuum limit at fixed m^L. 
This does not appear in our analysis, since we keep a physical mass fixed, when we take 
the continuum limit. 

The only remaining linear lattice artefact is cancelled by chosing 

ci' 1] = 0.0191410(1). (4.5) 

We conclude that the observable g~ 2 is 0(a) improved to one-loop order with our 
action (section 2.6) and the choices c sw = 1 and c t as quoted in eqs. (4.1), (4. 2) and 
(4.5). Furthermore, a similar analysis for v ^ shows that the same is true for the 
obervable v, eq. (2.51). At present, we do not know, however, whether Sf is the general 
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tree-level improved fermion action, since the effect of two or more fermionic surface 
operators could in principle cancel each other in our observables. The general structure 
of the 0(a) improved action deserves further investigations. 

4.3 Nonlinear lattice artefacts 

We now discuss the complete one-loop lattice artefacts of the step scaling function. 
Besides the case c, w = 1 and c[ 1A) = 0.019141, we also consider the situation without 
improvement (c sw = and cj 1 ' 1 ^ = 0). The latter may be of interest in Monte Carlo 
calculations because the Sheikholeslami-Wohlert term represents a significant numerical 
overhead in any simulation algorithm. 



L/a 


c sw = 1 and cj 1 ' 1 ^ = 
z = z = 1 


0.019141 

z = 2 


c sw = and 
z = 


.C 1 ' 1 ) n 
c* =0 

z = 2 


4 


0.00009 


0.00116 


0.00197 


-0.00273 


0.00035 


5 


-0.00005 


0.00069 


0.00126 


-0.00330 


0.00001 


6 


-0.00010 


0.00045 


0.00088 


-0.00346 


-0.00020 


7 


-0.00010 


0.00032 


0.00066 


-0.00344 


-0.00032 


8 


-0.00008 


0.00025 


0.00051 


-0.00334 


-0.00041 


9 


-0.00007 


0.00020 


0.00041 


-0.00322 


-0.00047 


10 


-0.00006 


0.00016 


0.00033 


-0.00309 


-0.00051 


11 


-0.00005 


0.00013 


0.00028 


-0.00296 


-0.00054 


12 


-0.00004 


0.00011 


0.00023 


-0.00284 


-0.00057 


13 


-0.00003 


0.00009 


0.00020 


-0.00273 


-0.00059 


14 


-0.00003 


0.00008 


0.00017 


-0.00263 


-0.00060 


15 


-0.00003 


0.00007 


0.00015 


-0.00253 


-0.00061 


16 


-0.00002 


0.00006 


0.00013 


-0.00244 


-0.00061 



Table 5: Lattice artefacts <5i,i in the step scaling function for various L/a. 



The lattice version of the step scaling function, defined exactly as in the continuum 
(cf. section 2.5) but for finite a/L, is denoted by S(s, u, z, a/L). Its lattice artefacts 
have the perturbative expansion 

S(2,u,z,q/X)- a(2,u,z) 2 

r = S 1 (n f ,z,a/L)u+0(u), (4.6) 

(T^Z, U, Z) 

with the one-loop coefficient 

t>i(nf, z, a/L) = 8 lt0 (a/L) + rifS ltl (z, a/L). (4.7) 

For the one-loop O(a) improved action, the lattice artefacts in the pure gauge theory 
vary from S lfi (l/6) = -0.004 to ^, (1/16) = -0.0004 (cf. table 1 in ref. [3]). Table 5 
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shows that even several massless flavors of quarks induce negligible lattice artefacts on 
this scale, when full 0(a) perturbative improvement is switched on. 

In contrast, each massless flavor introduces 8i t i(Q,a/L) ~ —0.003 lattice artefacts 
when we consider the action without improvement terms. Although the magnitude is 
not very large, there is a distinct problem. £1,1 varies little with the lattice spacing in 
the accessible range of a/L. Only around L/a = 16 does the linear decay of the lattice 
artefacts set in, a behavior which is dominantly caused by the coefficient si ^ in 
eq. (3.9). Clearly, lattice artefacts of this type are difficult to detect in a simulation. 

£1,1 behaves similarly for values of z up to z = 2. In particular, it decays roughly 
proportional to (a/L) 2 , when improvement is switched on. 

Very heavy flavors decouple in the continuum step scaling function up to 1/z- 
corrections. When one stays at fixed I, however, and increases z, one enters the regime 
arriQ > 1. It is evident that one needs to avoid this region, where all fermionic effects 
are accompanied by large cutoff effects. One might think, however, that - as in the 
continuum - such a fermion decouples and it does not really matter whether one keeps 
it in the simulation or not. We want to briefly explain that there is a subtlety, here. 

It is easy to expand pi^ for large amo (keeping I fixed) in order to see qualitatively 
what happens in this limit. As long as cj 1 ' 1 ^ = 0, the fermion contribution to £ starts 
with a term oc 1/(1 + amo) = exp( — z/l). It decouples faster than in the continuum. 
When we switch on O(a)-improvement, however, a term cj 1 ' 1 ^// remains, since a fermion 
with amo > 1 does not at all contribute the lattice artefact calculated in an expansion 
in a. Thus O(a)-improvement is ruined when amo > 1- Note that this behavior is not 
seen in table 5, since we restricted ourselves to amo < 1/2. 

Let us finally point out that the exact numbers listed in table 5 depend on the 
convention that we chose for the renormalized quark mass, eq. (3.4). Nevertheless, we 
do not expect the qualitative behavior of £1,1 to be very sensitive to the choice of m. In 
particular, the large effect of tree-level improvement is caused by the vanishing of si in 
eq. (3.9) with tree-level improvement. This is the case irrespective of the definition of 
m. 

5 Summary and conclusions 

This article presents one step in the effort to compute the strong coupling from first 
principles following the original idea of refs. [1,8]. The challenge is to compute non- 
perturbatively the evolution of a renormalized coupling from low to high energies. As 
such a computation is to be performed through numerical simulations on finite systems, 
it is non-trivial to cover a large range of scales and take the continuum limit at the same 
time. If the coupling is defined in finite volume, running with the linear dimension, L, 
of the system, this potential problem can be solved by stepping up the energy ladder 
recursively. A suitable definition of a coupling has been known for the pure gauge 
theory [8,3]. 
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Based on recent work on the QCD Schrodinger functional [9,10], we have extended 
the definition of the coupling to include quark fields in the functional integral. In par- 
ticular we have made use of the freedom that exists in chosing the spatial boundary 
conditions for the quark fields. Requiring periodicity up to a phase, 9, (cf. eq. (2.6)), 
this phase can be tuned such that the Dirac operator taken in the background field has 
a relatively large lowest eigenvalue. This occurs around 9 = x/5. In a numerical simu- 
lation in a sufficiently small volume, the gauge field fluctuates around the background 
field. Thus we expect that numerical simulations are eased by the corresponding choice 
of 9 = x/5. 

We computed the renormalized coupling at one-loop order using both dimensional 
regularization and lattice perturbation theory. Our main result is the fermion contribu- 
tion to the matching coefficient between the coupling in the Schrodinger functional and 
the MS scheme. It is listed in eqs.(3.15-3.18) and table 4. We found that the quarks 
generally introduce rather small effects. 

The Schrodinger functional provides us automatically with a mass-dependent cou- 
pling. When one crosses the threshold around i _1 = m, the /3-function changes rather 
slowly from its value for massless quarks to the large mass limit corresponding to the 
theory with that quark removed. One might be worried that such a slow threshold be- 
havior is difficult to resolve in a Monte Carlo simulation, because quarks of a large mass 
introduce significant cutoff effects. Investigating this question quantitatively, we found, 
however, that possibly resulting uncertainties are small on the level of the statistical 
errors of corresponding pure gauge theory results. Thus there is no reason to expect 
that the threshold behavior seen at one-loop precision poses a problem in the numerical 
simulations. 

As the final goal is to compute the evolution of the coupling in the continuum limit, 
it is important to investigate the size of lattice artefacts. They have been computed 
to one-loop accuracy for the step scaling function (2.35). We have considered two 
different lattice actions. I) The Wilson action, with the boundary terms taken exactly 
as in ref. [9]. II) The Wilson action plus the Sheikholeslami Wohlert improvement term 
removing 0(a) "volume-effects" and the surface term that compensates for the operator 
J2k=i tr{-Fbfe-fbfe} m Symanzik's local effective Lagrangian. The one-loop contribution to 
the coefficient of the latter is determined through the evaluation of the lattice artefacts 
that are linear in a. Its value is given in eq. (4.5). 

The results obtained with the two actions are as follows. The overall size of lattice 
artefacts introduced by one quark flavor into the step scaling function is rather small, 
namely they are below the per-cent level for couplings g 2 < 3 and reasonable lattices 
(L/a > 4). For the improved action (II) and massless quarks, they are even more than 
another order of magnitude smaller. Presumably the latter finding is rather accidental 
for our specific observable. A more important difference is the a/i-dependence of the 
lattice artefact for the Wilson action (I): at zero quark mass, it practically fakes a 
constant for all lattice sizes that are accessible to simulations. The slow vanishing of 
the lattice artefacts sets in only around L/a = 16. 
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Despite the fact that our explicit one-loop computations show that the step scaling 
function can be O(a)-improved by including the two operators mentioned above, the 
complete structure of the O(a)-improved action of the Schrodinger functional is not 
known, yet. It is possible that additional operators are necessary beyond one loop. 
Furthermore, even at the one-loop level we can not exclude that the improvement of 
other observables necessitates further operators in the action. This problem is currently 
investigated. 

A further important question is a non-perturbative definition of the renormalized 
quark mass that can be well computed in the numerical simulations and shows small 
lattice artefacts. Also this problem is presently being addressed. 

We would like to thank M. Liischer, P. Weisz and U. Wolff for helpful discussions 
and critical comments at various stages of this work. One of the authors (S.S.) acknowl- 
edges financial support by CERN, where part of this work has been done. 

Appendix A 

In this appendix we provide the details of our computation of pi^z, I) which makes use 
of a first order recursion relation for the lattice Dirac operator (cf. section 3). To increase 
the readability, we use lattice units (a = 1, and thus L = I) and assume summation 
over repeated indices throughout the appendix. 

The lattice operator D$ = ~f5(D + mo) (cf. eq. 2.40) is hermitian with real eigen- 
values that come in pairs ±\/An, with A n , (n = 0, 1, . . . , n m ax) being the eigenvalues of 
A2 (cf. table 2 for the numerical values of the first few A n ). Therefore, the determinant 
of D$ is positive, and the quark field contribution pi t i(z,l) to the one-loop coefficient 
Pi (3.1) can be written 



provided we have set c t = 1. Due to the special properties of the background field (cf. sec- 
tion 2), the eigenfunctions of the lattice operator D$ are of the form 



Pi,i{ z -> = (k n f) lndetZ>5 



j7 =l /=0' 



(A.l) 



ip(x) = exp(ip k x k )u nc f(x ). 



(A.2) 



Here, {u nc , n c = 1, 2, 3} denotes the canonical basis in color space, 




(A.3) 



and pk are the spatial components of the allowed momenta 



p k = (2xn fe + 0)1 L, As = 1,2,3, 



(A.4) 
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with integer numbers n k = 1, . . . , L. 

According to the structure of the eigenfunctions, the determinant of D$ factorizes 
and eq.(A.l) becomes 

fi,i(^o = ^EEi lndet ^U=o- ( A - 5 ) 

n c =l p ' 

The reduced operator V$ acts on the functions f(xo) (A. 2) which live in the subspace 
of fixed spatial momentum p and color n c . Setting xq = t, we have 

(V 5 f)(t) = - l5 P-f(t + 1) + l5 h(t)f(t) - l5 P + f(t - 1), (A.6) 

with the boundary conditions (2.5) 

P+/(0) = 0, P_/(i) = 0. (A.7) 

To write down the explicit expression for the coefficient function h(t), we introduce the 
notation 

q k (t) = u)t + r k , q k (t) = smq k (t), q k (t) = 2 sin \q k (t), A: = 1,2, 3, 

(A.8) 

where u) and r k are related to the color component n c of the boundary fields (2.9), 

u = (<j>' nc - K C )IL\ r k = Pk + <j> nc /L. (A.9) 

We then have 

1 3 1 
h(t) = 1 + m + - E ik{tf + iq k (t)j k - -c sw j j k p ok , (A.IO) 

1 k=i 1 

where pofe denotes the color component n c of the lattice field tensor (2.47) and further 
evaluates to pofe = * sin a;, independently of k = 1, 2, 3. 

At first sight, V$ seems to act like a second order difference operator. In fact, 
the projectors P± hide the first order structure which can be recovered through the 
following re-formulation. For < t < L we define a new function F(t), 

F(t) = P-f(t) + P+f(t-l), (A.ll) 

for which the boundary conditions (A.7) take the form 

P+F(l) = 0, P_F(L) = 0. (A.12) 

It is then straightforward to show that the eigenvalue equation, (D5 — /x)/ = 0, is 
equivalent to the recursion relation 

F(t + 1) = A(t)F(t), (A.13) 
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with A(t) being explicitly given by 

A(t) = -ait)' 1 |p_ [n 2 - a(t) 2 + /i 75 (c fe (f )7fc - b k {t) lk + l) 

+ c fc (t)7*(&j(*)7j-l)] (A.14) 
+ P+ [b k (t) lk - /i 75 - l] J. 
Here, the coefficient functions are 



1 3 

a(t) = l + m + -^ftW 2 , 

Z fe=i 

6 fe (i) = ig fe (f) - -c sw p ok , (A.15) 

c k (t) = iq k (t) + ^c sw p fe, A; = 1,2,3. 

If we now prescribe a value for F(l), we may calculate F(L) for any L > 1 by solving 
the recursion relation (A. 13), 

F(L) = A(L - 1)A(L - 2) • • • A(2)A(1)F(1). (A.16) 

At this point it is convenient to introduce the 2x2 matrix M(/x) in the subspace defined 
by the projector P_, viz 

M(/x) = (B(L -l)---B(l))__, B(t) = a(t)A(t). (A.17) 

The boundary conditions (A. 12) for F then correspond to the requirement 

M( M )F(1)_ = 0, (A.18) 

with the non- vanishing 2-vector _F(1)_. To have a solution to eq.(A.18) the determinant 
of M(/x) must vanish. One then notices that detM(/x) is a polynomial in /x of degree 
4(i — 1), and, using the same arguments as in appendix C of ref. [8], one concludes that 
detM(/x) is proportional to the characteristic polynomial of V$. In fact, the correct 
normalization has been anticipated in eq. (A.17), 

det(£> 5 - /x) = det M(/x), (A.19) 

so that one finally obtains 

— lndetPg = TrfM" 1 — M), M = M(0). (A.20) 

OT] OT] 

Note that the trace in the r.h.s. of this equation is over a 2 X 2 matrix. The inversion 
of the regular matrix M is trivial and its derivative with respect to tj can be calculated 
using Leibniz' rule for the product of matrices in eq. (A.17). 
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